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Abstract 



We consider the most general solar model, using the neutrino fluxes as free 
parameters constrained only by the solar luminosity, and show that the com- 
bined solar neutrino data exclude any astrophysical solution at 98% C.L. Our 
best fit to the ''Be and fluxes is respectively <7% and 37ib4% of the stan- 
dard solar model prediction, but only with a large (5-6 for 1 d.f.). This best 
fit to the fluxes contradicts explicit nonstandard solar models, which generally 
reduce the ^B flux more than the ^Be. Those models are well parameterized 
by a single parameter, the central temperature. 
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Each of the solar neutrino experiments show deficits of the solar neutrino flux 

compared to the standard solar model (SSM) predictions as summarized in Table I. 
Numerous astrophysical solutions have been proposed to explain the discrepancy between 
theory and experiments p|-p!2|. One category of such proposals changes the input parameters 



of the solar models, assuming that the uncertainties of those quantities might be significantly 
underestimated in the SSMs. For example, the neutrino fluxes are known to be sensitive to 
the opacity, and explicit models have been constructed with smaller values of the opacity or 
smaller values of the heavy element abundance. The nuclear reaction cross sections, which 
are extrapolated from laboratory conditions, are another potential source of uncertainties: 
there might be some mechanism affecting the low energy cross sections and reducing the 
neutrino production, and such effects might even be correlated among the different reac- 
tions. A second category of proposals attributes the neutrino deficit to mechanisms such 
as rotation, magnetic fields, turbulent diffusion, mixing of elements, or hypothetical weakly 
interacting particles (WIMPs) that are not included in the SSMs. Both kinds of theoretical 
proposals usually reduce the expected flux of high- and medium-energy neutrino production 
by lowering the temperature in the core region where the nuclear fusion takes place, and 
can be parameterized by a lower central temperature (T^). In previous studies using the 
power law dependence of the neutrino fluxes on Tc [[l3i |lH3, it was shown quantitatively 
that such cooler sun models are incompatible with the experimental data, especially because 
the higher Kamiokande observed rate relative to the Homestake rate cannot be explained so 
long as cool sun models reduce the expected *B flux more than the ''Be flux. This failure of 
astrophysical resolutions of the solar neutrino deficit suggests particle physics solutions, such 
as the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism [|15|, which fits all observations 
and is taken as a strong hint of neutrino mass and mixings |I^p!6| -pO|] . 

In this paper, we remove the assumption of a power law dependence and examine ar- 
bitrary solar models by allowing the four relevant neutrino fluxes (f){pp), 0(Be), 0(B), and 
0(CNO) to change freely. We do not advocate such models or claim that they are consis- 
tent with other solar observations, but only show that they are incompatible with the solar 



neutrino data. In our most general solar model, we assume: (1) The Sun is in quasi-static 
equilibrium and generates energy by nuclear fusion in the pp and the CNO chains; (2) As- 
trophysical mechanisms may change the magnitude of each neutrino flux component, but 
do not significantly distort the energy spectra of the individual components (Particle 
physics solutions such as the MSW effect often depend on neutrino energy and therefore do 
distort the neutrino spectrum.); (3) All reported experimental results are correct, as well 
as the calculations of radio-chemical detector cross sections. Because the Kamiokande and 
Homestake results are crucial to our conclusions, we will also consider the possibilities that 
their uncertainties have been underestimated. 

By our first assumption, the well-measured solar luminosity imposes the constraint 

(p{pp) + 0.979 0(Be) + 0.955 0(CNO) = 6.57 x 10^° cm'^s-^ (1) 

among the pp, ^Be, and CNO fluxes, when the different energies carried off by neutrinos are 
taken into account. 

In Figs. we present the results of all solar neutrino experiments in the 0(Be) — 0(B) 
plane in units of the Bahcall-Pinsonneault predicted fluxes. Essentially all astrophysical so- 
lutions, including insensible models, are represented in the plane, from (f)(Be)/(f)(Be)ssM = 
0(B)/0(B)5SAf = 1 for the SSM to 0(Be)/0(Be)s5M = 0(B)/0(B)s5m = for the minimum 
rate model p2| . Fig. |I] shows the constraints from each experiment obtained by minimizing 



the with respect to (p{pp) and 0(CNO) at each point subject to the luminosity constraint. 
Our fit includes experimental uncertainties as well as detector cross section uncertainties. 
The uncertainties of minor fluxes {pep, hep, and ^^F) are included, but contribute negligibly. 
The Kamiokande result constrains only the ^B flux, while the Homestake data and combined 
SAGE-GALLEX data constrain both the ^Be and ^B fluxes. At 90% confidence level (C.L.), 
none of the experiments are consistent with the Bahcall-Pinsonneault SSM. Indeed, the com- 
bined experiments together allow only a small parameter space around (l)(Be)/(j){Be)ssM ~ 
and (f)(B)/(f)(B)ssM ~ 0.4, but with a large value: 5.6 for degrees of freedom (3 data 
- (4 parameters - 1 constraint)). No general statistical interpretation exists in such a case. 



other than to conclude that this model is excluded. If one considers the Be flux to be fixed 
at 0, then the fit has 1 degree of freedom and this possibility is excluded at the 98% C.L. 
This shows that any astrophysical solution in which the spectral shape of the individual 
neutrino fluxes is unchanged is incompatible with observations. 

Fig. shows the confidence levels of the combined fit in the two dimensional 0(Be) —0(B) 
subspace. The contours are determined by = Xmin~^ '^X^ : where Xmin obtained allowing 
an unphysical negative ^Be flux. We could alternately have taken Xmin physical region 

by restricting the probability distribution to the physical region (0(Be) > 0). This procedure 
would have ignored the fact that the best fit is very poor, and would grossly overestimates 
the allowed region. We therefore present our results as a qualitative display of the confidence 
levels. 

We also show in Fig. |^ the Bahcall-Pinsonneault SSM with 90% C.L. uncertainties, the 
1000 Monte-Carlo SSMs of Bahcall and Ulrich [Q, the central value of the Turck-Chieze- 
Lopes (TC) SSM 0, and various explicit nonstandard solar models constructed to solve the 
solar neutrino problem: the low Z model in which the heavy element abundance is reduced 
by 90% from the standard value the low opacity models with 10 and 20% reduced opacity 
0; the solar models with increased pp cross sections (^n) by 30, 50, 80, 100, and 150% 
from the SSM value [|l^; and the solar model with WIMPs ||T2|,|8|. Also the power laws for 
the core temperature and 6*11 obtained from the Monte-Carlo SSMs are extrapolated from 
the SSM region and displayed. The uncertainty due to the p+'^Be cross section (9.3%) is 
shown as error bars. 

Because the decay of ^B follows the reaction p +^ Be — >^ B + 7, all explicit nonstandard 
models predict more reduction of the ®B flux than the ^Be flux (i.e., 0(Be)/0(Be)s'5M > 
(f){B)/(f){B)ssM)- Any reduction of the ^Be production rate affects both the ^B and ^Be 
flux equally. Other uncertainties in the p+^Be rate affect only the ^B flux. Therefore, 
unless there is some independent mechanism to suppress only the ''Be neutrino emission, 
all realistic nonstandard solar models are in serious contradiction to the solar neutrino 
data, which constrain the two fluxes to (f){Be)/(j)(Be)ssM < 0.07 and (j){B)/(f){B)ssM = 
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0.37 ± 0.04 (Icr). This emphasizes that there are two solar neutrino problems: (I) The 
neutrino fluxes observed in every experiment are significantly below SSM predictions at 
90% CL. (11) Kamiokande and Homestake together allow only the very implausible fit 

0(Be)/0(Be)55A/ « 0(B)/0(B)55m. 

The two curves in Figs. assume that the ''Be and ^B neutrino fluxes each depend 
simply on powers of the central temperature Tc (solid curve) or of the pp nuclear cross section 
factor 5*11 (dot-dashed curve), while the pp and CNO neutrino fluxes are adjusted to obey the 
solar luminosity constraint (Eq. 0). For the solid curve we assumed 0(Be) ~ T®, 0(B) ~ T^^ 
so that 0(B) oc 0(Be)^'^^. For the dot-dashed curve we assumed 0(Be) ~ Sf^^'^^, 0(B) ~ 
Sfi^-^^ so that 0(B) oc 0(Be)2-67. Those exponents were obtained by Bahcall and Ulrich from 
1000 SSMs with input parameters randomly distributed near the most probable values ||^. 

The nonstandard solar models (the low opacity models, the low Z model, and the models 
with large Su) illustrated in Fig. ^ include physically unreasonable models with Sn as 
large as 2.5 times and Tq as small as 0.97 times their most probable values. Within their 
theoretical uncertainties, all model predict 0(B) oc 0(Be)" with n = 2.25 — 2.67. Most ex- 
tremely nonstandard solar models still lead to ^Be and ^B neutrino fluxes that are adequately 
parameterized by simple power laws, e.g., in the central temperature or Sn. (Exceptions 
include the maximum rate model P], the WIMP model, and the model with Su = |^.) 
This happens because, although the Sun as whole is not self-homologous (polytropic) , over 
the range of temperatures and densities in the Sun's inner core (r < 0.2Rq), 91% of the 
neutrino and energy production derives from the single pp reaction, and all the nuclear re- 
actions and opacities in the present Sun can be approximated by power laws. Consequently, 
when the luminosity is held constant, large changes in input parameters lead only to nearly 
homologous changes in core temperatures, mass, and radius. 

So far, we have shown that the Kamiokande and Homestake results together, if correct, 
essentially exclude any astrophysical solutions. What if either experiment were wrong? In 
Fig. 1^, we show the enlarged allowed region of the combined flt when Homestake's quoted 
experimental error is tripled. The data still strongly disfavor the nonstandard solar models: 
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the cooler sun with Tc reduced by 5% is only allowed at ~1% C.L. Expressed otherwise, 
the best fit with 0(Be) = corresponds to a Homestake rate of 2.9 SNU; the best cool sun 
fit when the Homestake error is tripled is 3.3 SNU, compared with the value of 2.23 ± 0.23 
in Table I. We have also carried out a calculation with the cross section uncertainties for 
the chlorine and the gallium detectors increased by factors of three, and obtained a similar 
result (Fig. Of course, if one entirely disregards either of these two experiments, a large 
class of nonstandard models become possible. 

In summary we have considered the most general solar model with minimal constraints 
using the neutrino fluxes as free parameters, and shown that the fit is excluded by the 
solar neutrino data at 98% C.L., i.e., essentially any astrophysical solution is incompat- 
ible with the quoted data. Furthermore, this very improbable best fit point requires 
0(Be)/0(Be)55M < 0.07 and (j){B)/(t){B)ssM = 0.37 ± 0.04 (la), which is inconsistent with 
virtually all explicit nonstandard solar models, which predict a larger reduction of the 
flux than the ''Be flux. Increasing the Homestake experimental error or the detector cross 
section errors by factors three does not justify the nonstandard solar models. 

We conclude that at least one of our original assumptions are wrong, either (1) Some 
mechanism other than the pp and CNO chains generates the solar luminosity, or the Sun 
is not in quasi-static equilibrium; (2) The neutrino energy spectrum is distorted by some 
mechanism such as the MSW effect; (3) Either the Kamiokande or Homestake result is 
grossly wrong. 

We also noted that almost all explicit nonstandard models fall on a narrow band in 
the 0(Be) — 0(B) plane, and can be characterized by a single effective parameter, the core 
temperature [|13],|14|. 

It is pleasure to thank Eugene Beier for useful discussions. This work is supported by 
the Department of Energy Contract DE-AC02-76-ERO-3071. 
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TABLES 

TABLE L The standard solar model predictions of Bahcall and Pinsonneault Q and of 
Turck-Cliieze and Lopes 0, along with the results of the solar neutrino experiments. 

BP SSM TCL SSM Experiments 

Kamiokande 1 ± 0.14 0.77±0.19 0.50±0.07 BP-SSM 

Homestake (CI) 8±1 SNU 6.4±L4 SNU 2.23±0.23 SNU (0.28±0.03 BP-SSM) 

SAGE & GALLEX (Ga) 13L5l^ SNU 122.5±7 SNU 71±15 SNU (0.54±0.11 BP-SSM) 
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FIGURES 

FIG. 1. Each experiment is fit to the pp, "^Be, ^B, and CNO fluxes, imposing only the 
luminosity constraint. The fit neutrino fluxes are plotted in the "^Be-^B plane in units of the Bah- 
call-Pinsonneault predicted fluxes p. This parameter space represents all possible astrophysical 
solutions consistent with our (minimal) assumptions. The 90% C.L. uncertainties of the Bah- 
call-Pinsonneault SSM are shown in the upper-right corner. 

FIG. 2. The allowed region from the combined fit of the Kamiokande, Homestake, and gallium 
results at 90, 95, and 99% C.L. allowing (unphysical) negative values for 0(Be). For (j)(Be) > 0, 

= 5.6 for the best fit, and the model is excluded at 98% C.L. for 1 degree of freedom. Therefore 
any astrophysical solution is excluded at > 98% C.L. Also shown are various nonstandard solar 
models and the power laws of the central temperature (Tc) and the pp cross section (Su) obtained 
from the Bahcall-Ulrich SSMs and extrapolated from the SSM region. All nonstandard models 
other than WIMPs can be approximately parameterized by Tq or Su. The error bars show the 
uncertainty of (j)(B) due to the p+^Be cross section [^. 

FIG. 3. The combined fit when the Homestake experimental error is tripled. 

FIG. 4. The combined fit when the detector cross section uncertainties are tripled from 3.3% 
(Homestake) and 4% (gallium). 
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